Finite Size Effect in Persistence in Random Walk 



o 
o 

(N 
D 



D.Chakraborty * and J.K.Bhattacharjee^ 

Department of Theoretical Physics, 
Indian Association for the Cultivation of Science, 
Jadavpur, Calcutta 700 032, India. 

We have investigated the random walk problem in a finite system and studied the crossover 
induced in the persistence probability by the system size. Analytical and numerical work show that 
the scaling function is an exponentially decaying function. We consider two cases of trapping, one 
by a box of size L and the other by a harmonic trap. Our analytic calculations are supported 
by numerical works. We also present numerical results on the harmonically trapped randomly 
accelerated particle and the randomly accelerated particle with viscous drag. 



I 

o 
o 



(N 

> 

cn 
o 
o 



X3 
O 

o 



X 
J3 



The phenomenon of persistence has attracted a lot of 
interest in the recent years, both theoretically Jj-j)] as 
well as experimentallv[To | -[l ^ . The word "persistence" 
itself conveys the meaning of survival. Associated with 
this survival is the survival probability It is simply 
the probability that the local field has not yet changed 
its sign upto time t. For a wide range of models the 
survival probability decays as a power law, that is 
p{t) ~ , where is a new non-trivial exponent called 
the persistence exponent. Established results, exist 
for many models- Random Walk problem. Diffusion 
problem l!i,Ising Model with Glauber dynamics 0, 
Surface growth 6] , Phase ordering kinetics j5j . 

In an experimental setup, however, finite size effects ap- 
pear because of the size of the apparatus and boundary 
effects come into play in the dynamics of the system. As 
a result the survival probability also depends on the fi- 
nite size parameter. A particularly clear example of the 
crossover effects induced by finite size is the recent ex- 
periment on a single polystyrene sphere in a harmonic 
potential 0| It is this demonstration of crossover effect 
that has motivated us in our present work in investigat- 
ing the finite size effect on the survival probability and 
how it scales with the finite size parameter. 

There has been an investigation of this for the Ising 
system in higher dimension [20| which lead to the con- 
clusion that finite size scaling in the usual sense holds for 
persistence as well. The system that we investigate here 
do not show finite size scaling in exactly that sense, al- 
though they exhibit pronounced finite size effects. We 
have considered an analytically solvable model in the 
present work- the case of a brownian particle confined 
in a box and the case of a brownian particle trapped 
by a harmonic potential. We find that the survival 
probability P{t, L) does not have the usual scaling form 
P{t, I) ^ t^^ f {jj), with f{x) constant for x — > and 
f{x) ~ for X » I. Instead, we find that P{t,L) can 



be expressed as P[t,L) — t^^J{^) but f[x) expo- 
nentially as a; — > oo. This is consistent with the generic 
form anticipated by Redner0|. The exponent z is found 
to be 2.0 as in the case of Manoj and Ray [2(tI |. 

The simplest of all the models for which there exists 
an established result is the Random walk problem. A 
random walker obeys a differential equation of the form 
dx{t) 



dt 



— Tjit) where ri{t) is a white noise. To find the sur- 
vival probability we ask the question whether the quan- 
tity sgn[x{t)— < x{t) >] has changed its sign upto time 
t. The survival probability Po(T) in terms of the variable 
a = sgn{X{T)) can be found from A{T) =< a{0)a{T) > 
[3. In this case Pq{T) = {2/Tr)sin-^[exp{-XT)]. For a 
random walker A = 1/2 and the survival probability goes 
as Po{T) ~ exp{—T/2) 18] . Analytical and numerical 
result show that the probability goes as p{t) ~ t~2 and 
the persistence exponent in this case is 9 = ^ [l^- 

We have investigated the finite size effect in the 
random walk problem in two ways. Firstly, the random 
walker is constrained to move in a box with reflective 
boundaries at a; = ±L The probability distribution 
P{x, t) in this case obeys a diffusion equation with 
an appropriate boundary condition. Solution to the 
diffusion equation with the proper boundary condition 
gives P{x, t). In the second problem, the random walker 
is trapped in a harmonic potential. Both the problems 
are analogous to each other with the identification 
w ~ In both cases we calculate the correlator 

0(^17^2) x{t2)x{ti) >, where x{t) is the value of x 
at time t. To make it a Gaussian Stationary Process 

xU) 

(G.S.P) we transform x{t) io X ^ ^ ' 



and a 



> 



suitable transformation for the time variable from t to 
T. Thus the correlator 0(^2, ii) ^ f{\T2 - Ti\). From 
the correlator /(T) we get the survival probability p{t). 

We first consider a particle in one-dimension perform- 
ing random walk. The equation governing the dynamics 
of the particle is given by 
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where x{t) is the displacement of the particle and r]{t) is 
a random function whose moments are given by 





DS{t~t') 



(2a) 
(2b) 



where D is the diffusion coefficient. In the present 
problem we confine the motion of the particle within a 
cage with boundaries at a; = ±L. The boundary of the 
cage is reflective that is upon reaching the boundary the 
particle is reflected to the nearest lattice site. 

The probability P{x,t) that the coordinate is a; at a 
time t starting from x ~ a.t t ~ obeys the diffusion 
equation 



dPixA) d^P{x,t) 



dt 



D- 



(3) 



with proper boundary condition. Since the particle is 
reflected from the boundary, the particle current x = ±L 
must be zero. Thus at the boundary 



-D 



dP{x,t) 



dx 







(4) 



x=±L 



The solution to Eq.© with the boundary condition 
Eq.Q) can be written as 



mix 

P[x,t) = cos(— — ) e 



(5) 



The complete form of the probability P{x, t) taking into 
consideration the normalization can be written down as 



P{x,t) 



I 1 



COS 



n=l 



WKX 



(6) 



The average of the square of the displacement is given by 



< (Ax)2 > = — + 



4i2 ^ 1 
— cos[mT)e 



TT 



n=l 



(7) 



This is the exact answer. As expected it exhibits finite 
size scaling and can be cast in the form 



t 



(8) 



such that g{j-! , ^ 
and g{{ 



1 for 



(i.e infinite system size) 
oo that is the extreme case of 



t 
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finite system. 

To see this the sum in Eq.Q can be decomposed as 



""FT" 



n odd 



I 



5Z ^7,2 



(9) 



In the limit L — > oo, all the modes in the summation of 
Eq.Q must be considered. It is then easy to see that 



< (Aa;)2 >- 



for t oo and L finite 



(10) 



and taking the opposite limit we find that 

< {/Sxf >-> 2Dt for L -> oo and t finite (II) 

as expected. Keeping in mind our future need where 
f{x) ^ I as X ^ but decays very fast for x >> L, 
we will express Eq.Q as an approximate which is easy 
to handle. This is done using the Euler-Maclaurin sum 
formula for the two sums in Eq. @ . 

Keeping only e from among the different expo- 

nential decays and working to O(-p-) the coefficient of 

e the crossover function can be written as 



< (Ax)^ >= 



i2 
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(12) 



A more drastic approximation yields the expression 



< (Ax)2 >= 



d?x(t)j-2 



1_( l + (^2_g)_\ 



(13) 

The part in the bracket is the approximation for the func- 
tion g{j2)- For -p- — > oo, Ea. (|I3(l correctly reduces to 
while for -p- << I, we expand Ea. (|I3|l to obtain 



< {/\xf >= 2Dt 



1-7 



Dt 

Z2 



where 



7 = 71-2(1 ) 

' ^ 12^ 



(14) 



(15) 



The leading term in Ea. HI4|l is the correct limit for the 
unbounded system and the second term is the first cor- 
rection for finite L. 

We now proceed to calculate the correlator < 
x{t2)x{ti) > for the dynamics of Eq.(^ keeping in mind 
the approximations used in arriving at Ea. HI3l) for < 
(Aa;(t))2 >. Considering the equation for the proba- 
bility distribution (Eq.Q, we write down the expression 
for P{x2,t2', Xi,ti) the probability of finding a value X2 
at t = t2 if the value was xi at t = ti. We note the exact 
result 



[x{t2) - x{ti)]'^P{x2,t2;Xi,ti)dx2dxi 

^< {Ax{t2 - h)f > (16) 
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whence 



< x'^{t2) > + < x'^ih) > -2 < x{t2)x{ti) > 

= < iAxit2 - h) f > (17) 

and we can now use Ea. (|13|l to calculate a(ti,t2) =< 
x{t2)x{ti) >. To obtain a Gaussian Stationary Process, 

it is necessary to calculate the correlation of X — f''^^ 

and this leads to complicated looking expression. To ex- 
press the final answer in a particularly simple form, we 
use the regime << 1 and then exponentiate to find 



^ 0,04979 



< Xit2)X{ti 



- e 

t2 




.. P(t2-tl) 



(18) 



This exponentiation is predicted by the generic form an- 
ticipated by Redner. In the process of our calculation, 
we find the numerical prefactors which are not present 
in the general arguments of Redner. For comparison 
with the numerical simulations, these prefactors are es- 
sential. We now perform the transformation in time 
t^T ^lnt + ^. The correlator /(T2, Ti)in the trans- 
formed variables is 



FIG. 1: Semi-log plot of t'ip{t) vs j^. The straight line is the 
best fit for the linear part of the curve. 



Ref |l9| except there the inertial effect cannot be ignored. 
The equation of motion is 



dx{t) 2 



< X{T2)X{Ti) e-i/2(T2-Ti) 



dt 



-\- u> X = rj(t) 



(23) 



(19) 



The process is now a Gaussian Stationary Pro- 
cess (G.S.P). The survival probability is now given by 



where 77 (i) is the random noise whose moments are given 
by Ea. lPajl and Ea. (|2bjl . The expression for x{t) then 
becomes 



p(t)=t-'/' =t-'"f{ — ) (20) 

To test Eq. H20|l , we have calculated p{t) numerically. 



x{t) = e-"'* / e'^'''ri(t')dt' 



The correlator < x{ti)x{t2) > is given by 



(24) 



— '^^^ 

e 



(21) 
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We expect the semi-log plot of t ^^'^pit) vs — to be 



< x{h)x{t2) > = e-^ / / e- 

Jo Jo 

< V{t'i)v{t'2) > dt[dt'2 (25) 



straight line with a slope of 7 = 
obtained from our calculations is 



7/2. The value of 7 Using Ea. l|2b(l we have 



— (1 ) = 0.8761 

2 ^ 12^ 



(22) 



< x{h)x{t2) >= ir^ie'"" - e-'" (26) 



The correlator in the new scaled variable x{t) X{t) = 



Numerical simulation of the process was done using vari- 
ous values of L. The probability was obtained by averag- 
ing over 10^ configurations. Numerically obtained value 
of 7 is 0.9482. This discrepancy can be attributed to the 
approximate form of Ea. (|18|) . A semi- log plot of t^pit) 
vs is shown in Fig 1. This clearly shows the validity 
of Redner's generic form and the reasonableness of our 
approximations in arriving at the numerical value for 7 
and the fact that z = 2. 

We next consider a particle trapped in a harmonic po- 
tential and acted upon by a random force. Instead of 
sharp boundaries, we now have the potential confining 
the particle. If the confining length is L, then on dimen- 
sional grounds we expect uj ~ 1/L This is the set up of 



^(«) 



has the form 



sinh{bj'^t2) 
sinh{uj^ti) 



\/<^'(*)> 

<Xiti)X{t2) >=e-"^('i-*^) 

Writing = —^e^ sinh{uj^t) we have 

f{Ti,T2) =< X{Ti)X{T2) e-^(^^-^=) 



(27) 



(28) 



where A = i. The process is now a gaussian stationary 
process. The survival probability can now be written as 



piT) = e 



-AT 



(29) 
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FIG. 2: Log-log plot of p{t) vs lo t. The line shows the func- 
tion .'^ ^ 
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FIG. 3: log-log plot of survival probability vs time. The solid 
line is the best fit line. The estimated persistence exponent 
by best fit is 6^ = 0.2502. 



and in real time the survival probability is 
1 



P{t) = 



sinh{uj t)] 



H/2 



1 



where 



sinh{x) 



(30) 



(31) 



As stated before z — 2 and f{x) ^ 1 as a; ^ 0, while 
f{x) — > as X — > oo. For w — s- the above expression for 
probability reduces to the normal random walk problem 
and the probability p{t) goes as 2 . The numerical data 
is obtained for three values of w. For t < XjuP' the es- 
timated value of the exponent by fitting the log-log plot 
with a straight line is found to be = 0.5055. 

Finally we present numerical studies of the survival 
probability for a harmonically trapped randomly accel- 
erated particle and for randomly accelerated particle with 
viscous drag. The persistence exponent for the randomly 
accelerated particle is = 0.25 The equation of 

motion for the particle is 



+ LiP'X = 77(t) 



(32) 



where r\(i)) is a gaussian white noise with correlator given 
by Ea. (|2al) and Ea. (|2b|l . Rescaling t = ujt we see that 
for T<<lort<<l/w the first term domainates and 
the equation is motion is that of a randomly acceler- 
ated particle. A plot of survival probability vs time is 
shown in Fig. 3. With uj behaving noted above, 

this corresponds to z = 1 for the dynamic component. 
The survival probability is obtained by averaging over 
10^ configurations. 



For the randomly accelerated particle with viscous 
drag the equation of motion is 



(fix dx , , 



(33) 



This is the crossover which is of particular interest in 
the experiment of Ref . As can be seen from our re- 
sults the cross over in the form of the correlation function 
occurs at i ^ jT. In this case however two regimes ex- 
ists. For t « 1/r the equation of motion is that of a 
randomly accelerated particle with the first term dom- 
inating and for t » l/F the second term domainates 
and the equation of motion is that of a random walker. 
Thus the survival probability also shows a crossover from 
the randomly accelerated regime to random walk regime. 
The estimated values of the exponents in the two regimes 
is tabulated below. 



value of r 


e for i < 1/r 


61 for i > 1/r 


0.8 


0.2638 


0.4970 


0.1 


0.2510 


0.4888 


0.01 


0.2538 


0.4797 



Numerically obtained values of survival probability is 
plotted against time in Fig. 4. The survival probabilities 
is obtained by averaging over 10^ configurations. One of 
the important finding of Ref (1&] was that in the realistic 
situation of the experiment the crossover in < (Ax)^ > 
occurred for t » j;. This was attributed to a mem- 
ory dependent damping term. In a future work, we will 
explore the effect of it on the persistence problem. 
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